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The multifractal formalism characterizes the scaling properties of a physical density p as a function 
of the distance L. To each singularity a of the field is attributed a fractal dimension for its support 
/(a). An alternative representation has been proposed by Jensen considering the distribution 
... of distances associated to a fixed mass. Computing these spectra for a multifractal Cantor set, it 

CO ' is shown that these two approaches are dual to each other, and that both spectra as well as the 

' moment scaling exponents are simply related. We apply the same inversion formalism to exponents 

, obtained for turbulent statistics in the GOY shell model and observe that the same duality relation 

(N . holds here. 
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Initially motivated by the statistical characterization of velocity fluctuations in turbulence ^ , the multifractal 0, Q 
formalism has been shown to be a powerful way of analyzing a large body of different problems. It provides a simple 

Uand elegant way of performing a "dimensional analysis" of singular fields. In turbulence, this approach has been 
applied to the fluctuations of the velocity field, and deviation from the sirnple (monofractal) Kolmogorov "El scaling 
^ ' of moments of different orders has been observed both experimentally 00,1^, 9J and numerically 10, 11] . It has been 
used also to characterize the growth probabilities of Diffusion-Limited Aggregation |l3 . Random resistor networks 
[T^ at the onset of percolation have also been studied using this formalism. Extension to damage and fracture 
models have been proposed. 

In those examples a local physical quantity m — referred to as a "mass" in the following for concreteness — is 
^ \ distributed in space (or time) and the formalism allows to characterize the statistical distribution of this quantity, or 
0^ ■ equivalently its moment of any order, as function of the system size L (or time interval) over which it is considered. 
The field is decomposed, according to its singularities, a, into a continuous set of fractal supports. The corresponding 
fractal dimension /(a) as a function of the singularity a of the field is the multifractal spectrum. Hence, the number 
' n(m) of elements of mass m such that 

y . rn ~ (1) 

scales as 

C ■ n{m) ~ (2) 

> 

From this function, the scaling of any statistical moment of the field can be computed. Defining the moment of 
order g, Mg, and its scaling with the system size as 



M, = ^mf ^L-(<^) (3) 

i 

we can relate the scaling exponents r(g) to the multifractal spectrum, through a simple Legendre transform |^ 

9 = -/'(«) (4) 
T{q) = + /(a) 

Recently, Jensenflj proposed to consider an alternative approach to characterize the same fields. Instead of studying 
the statistical distribution of mass m{L) over a fixed distance L, he proposed to consider the distribution of distances 
L{m) such that a fixed mass m is contained in each subset. From the initial description (m(L)'^) cx i"^'*) a naive 
expectation would have been that {L{mY{q)) oc m'. However, considering the GOY model[i^[l3| (as a toy-model 
for turbulence), it was shown that the latter expectation was violated 0. Instead a different scaling was observed 

{L[mY (X m (5) 
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FIG. 1: Construction of the multifractal Cantor set. Starting from an interval of length L and mass M (generation 0), the first 
generation is obtained by splitting the interval into three equal length segments, and dropping the middle part. The mass M 
is distributed into two unequal parts rM and sM respectively for the left and right interval. The same procedure is repeated 
recursively onto each interval. 

But apparently, the series of exponent 9{p) seemed unrelated to the T{q). This unexpected feature suggested to use 
this new scaling as a complementary statistical property of turbulence. 

In the following we will consider a simple example of a multifractal set, using the standard Cantor set construction, 
but endowing each interval with a different mass [3. This simple case study allows to obtain a direct evaluation 
of the two multifractal spectra, as well as the corresponding scaling exponents. We show that within this example 
both approaches are related through simple duality relations. We then discuss the applicability of the previously 
derived duality relations to the case of turbulence. Numerical estimates of the scaling exponents of length moments 
for fixed velocity differences are obtained from She and Levesque proposed form for the velocity moments scaling 
exponents. We can apply this formula to obtain the series of exponents for the inverse statistics. To compare these 
"static" data to more realistic dynamical turbulence data we extract the scaling exponents for forward and inverse 
statistics of the GOY shell model. This was already done in P but here we extend the analysis to negative values of 
the moments of the standard forward structure functions. We then apply the inversion formula and compare to the 
exponents obtained by direct measurements of the inverse structure function. We obtain quite good results in the 
comparisons of these data sets as will be discussed in particular in Section VI. Although it is by no means a proof, 
it gives an indication that at least in some cases, the inversion formula we derive (which has been previously been 
derived in other contexts as we discuss) gives a relation between the exponents of forward and inverse statistics. 

II. STANDARD MULTISCALING FOR THE CANTOR SET 

The interval of length L is split in three equal parts and the middle one is removed. The mass M is split in two 
unequal parts, rM and sM such that r + s ~ 1. After N repetitions of this procedure, we obtain a generation N 
structure. The size of each piece is £ = 3^^L. Its mass is m = r'^s^~^M, where i is the number of r choices leading 
to a specific part. The number of such intervals carrying the same mass is n = while the total number of parts 



is5 = E. (;) = 2^ 



N' . . . 

We go to the continuum limit and define the real x as i — xN. Using Stirling formula we have 



= 2 



~N N' 



= {x-{l-xY'--Y 

In order to bridge this computation with the standard way of defining the multifractal spectrum Q , we introduce 

_ log(m/A/) 

log(£/L) /^N 

jy^) - ~iog(e/L) 
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FIG. 2: Multifractal spectra for the Cantor set. The bold curve is the f{a) function, while the dotted curve shows the dual 
spectrum g{a). 



and we define ao = — log(s)/ log(3) and qi = — log(r)/ log(3). A simple algebra leads to 

(ai -ao)log(ai -ao) - (ai -Q;)log(Q;i -a) - (a-Q!o)log(Q!-ao) 
^^"^ = (ai-ao)log(3) 



III. MOMENT SCALING 



The moment of order r of the mass distribution is defined as Aq{N) = n{i){m{i)/M)'^ . It obeys the recursion 
formula 

Ag{N) = {r" + s'^)Ag{N - 1) (9) 
For the definition of the scahng exponent T{q) as Ag{N) oc {£/L)~^'-'^^ we can write 

^(')-^ 

One basic property of the muhifractal formahsm is that the scahng exponents T(q) can be related to the multifractal 
spectrum through a Legendre transform. Indeed, the moment can be evaluated as Ag{N) = Ea(^/-^)^^"^"^'"- Hence, 
T{q) = maxa[/(a) + qa\. This defines the strength of the singularity a which contributes dominantly to the moment 
of order q. 

q = -f'{a) 

T{q) = f{a) + qa ^ ' 

The symmetry property of the Legendre transform allows to express the reverse transformation as 

I f r ^ ^ ^'/f (12) 

1 fW = nq) -qa ^ ^ 



IV. CONSTANT MASS ENSEMBLE 



Now we introduce the alternative approach of truncating the hierarchical construction at a fixed mass and not fixed 
generation (or length). The mass m/M is chosen, and thus i and j are related by 

r^fii = m/M (13) 
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FIG. 3: Scaling exponents r(g) for the mass as a function of the interval length in bold line, and the dual scaling exponent 
9{q) in dotted line. 



The length of such intervals is 



- = 3-('+-'') (14) 
L 



and their number is n = ( . * ) . 

To mimick a similar construction as previously we define £/L = [m/M)^ , and n = (m/M)^^^-*. A simple computa- 
tion leads to 

_ /3(ai - gp) log(/3) -f /3(ai - ap) log(ai - gp) - (1 - /3ao) log(l - /3ao) - (/3q;i - 1) \og{(3ai - 1) 

(ai - Q;p)log(3) 

Comparison with the original multifractal spectrum shows that they are related through 

5(/3) = Pf{l/P) (16) 
In fact this key relation can be simply derived by noting that 

n =(£/L)/(") 

= (m/Af)/3/(i//3) (17) 
= (TO/M)ff('3) 

from which Eg 1161 results. 



V. MOMENT SCALING 



We introduce similar moments in the dual ensemble Bp — L){£/ LY and define their scaling exponents 9{p) 

as Bp oc {m/My^P\ As for the primal ensemble, the scaling exponents 9{p) are related to the multifractal spectrum 
g{P) through a Legendre transform 



and reciprocally 



P =-9'{P) /-1 8^ 



g{f3) ^e{p)-p(3 ^'^^ 
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Let us now use the duality relation^] to relate the two series of scaling exponents. 



P = -g'm 

= -/(«) + /'(«)« (20) 
= —T[q) + ra — qa 
= -r{q) 



and 



d{p) =9if3)+pP 

= {f{a)+p)/a ^21) 
= iT{q) - qa - T{q))/a 

= -q 

Thus the two functions t and 6 are related by a mere inversion up to sign reversals (see also |27j'). 

-e{-Tiq))^q (22) 

One can also visualize the above result by noting that the graph of 6{p) is obtained from that of r(q) through a simple 
symmetry with respect to the line passing through the origin and of direction (—1, 1). 

Therefore, contrary to what was initially proposed, the two series of exponents are not independent. They are 
linked by a duality relation, as the two multifractal spectra / and g. 



VI. APPLICATION TO TURBULENCE STATISTICS 



In order to test the relevance of the above analysis to a more physical application than the Cantor set, we resort to 
the framework of turbulence which was the initial context of the suggestion of these dual quantities. In this context, 
the local physical quantity of interest is the velocity fluctuation Au (instead of the mass in the above example) studied 
over a distance r which plays the role of the local scale £. The early suggestion by Kolmogorov of the scaling moments 
T{q) — q/3 was shown to break down due to intermittency thus defining a non-trivial series T(g) which has resisted all 
theoretical attempts to compute them up to now. Nevertheless, semi-empirical formulas have been proposed which 
accounts rather precisely for the numerical values of these exponents as determined experimentally. In particular, the 
She and Levesque (LlQj formula appears as an accurate fit. They proposed 

T{q) =q/9 + 2 {l-{2/3r/^) (23) 

Table gives the corresponding r exponents for selected positive moments. In order to extract the corresponding 
9 exponents for the dual statistics, we apply the inversion formula H22|) . Indeed, we then need to employ the She- 
Levesque formula for negative moments. From the derivation of this formula, it is not obvious that such extension is 
allowed by the assumptions made by She and Levesque. Nevertheless, we take the liberty to continue the formula to 
negative moments and obtain the list of exponents 6 listed in Table I. 

In order to compare these results with data obtained from numerical simulations of turbulence models, and even 
more importantly to test the inversion formula (|22|l , we turn to dynamical turbulence generations by shell models |10| . 
This kind of data were already applied in the paper by one of us which proposed the inverse structure functions . In 
turbulence theory, it is well know that scaling behavior of velocity field u(x, t) and the understanding of intermittency 
effects in fully developed turbulence is described in terms of standard structure functions defined as 

{Au^iey) ~ (24) 

where the difference is 



Au^{£) = u(x + r)-u(x) , £ = \r\ (25) 

The average in Ea. (|24|l is over space and time. We have assumed full isotropy of the velocity field. The set of 
exponents Tq forms a multiscaling spectrum jl9f . 

The corresponding dual structure functions is defined by considering the following quantities 

{l{Au^r) ^ \Au^f^ (26) 
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FIG. 4: A plot of the ordinary structure function from shell model data with moment q = —2, i.e. ((Amx(^))') versus the 
length scale £. Note the three different regimes: The small scale regime of the smooth behavior; the inertial scaling regime and 
outer cut-off. 



where the difference Amx is again defined as in Eq. H25|l and .^(Aux) is understood as the minimal distance in r, 
measured from x, for which the velocity difference exceeds the value Aux lij]. In other words, we fix a certain set 
of values of the velocity difference Amx. Starting out from the point x, we monitor the distances £{Au^) where 
the velocity differences are equal to the prescribed values. Performing an average over space and time the inverted 
structure functions Eq. (|26|l are obtained. 

The turbulence data are obtained from simulations of the GOY shell model [lol IT^ [l7| . This model is a rough 
approximation to the Navier-Stokes equations and is formulated on a discrete set of fc- values, fc„ — r". We use the 
standard value r — 2. In term of a complex Fourier mode, u„, of the velocity field the model reads 

^<-l<-2) + fSnA, (27) 

with boundary conditions bi — bpf — ci ~ C2 ~ ajv-i = o,N = 0. / is an external, constant forcing, here on the forth 
mode. The coefficients of the non-linear terms must follow the relation a„ -I- bn+i + c„+2 = in order to satisfy the 
conservation of energy, E = |u„|^, when / = z/ = 0. The constraints still leave a free parameter e so that one 
can set a„ = 1, 5„+i — ~e, c„+2 = — (1 — e) [IDl- observed in one obtains the canonical value e — 1/2, if 
helicity conservation is also demanded. The set H27I) of N coupled ordinary differential equations can be numerically 
integrated by standard techniques. We have used standard parameters in this paper N — 21,v = 10~^,fco — 0.05, 
and / = 5 • 10"^. 

The structure functions exponents Tq are shown in Fig. Elfor integer moments in the interval q G [—10; 12]. A line 
connects the points in order to guide the eye. The associated exponents 9{q) for the inverse structure functions are also 
shown in Fig. [Slfor moments in the interval q G [0; 12], and are connected by a dotted line. It is possible to extract the 
exponents Tq to reasonable accuracy for negative moments —q although the quality of the scaling gradually decreases 
with the value of q. As an example we show in Fig. 2|thc behavior of (£(Aux)~^), corresponding to q — —2. We 
observe three distinct regimes: the small scales referring to the trivial smooth regime, the "inertial" scaling regime and 
the cut-off regime at large scales. For increasing value of negative moments, the point where (£(Aux)'') are small will 
be enhanced. This is a very important point to be considered for the analysis of experimental data. Indeed, the lower 
cut-off of the inertial, increasing with large negative values of q 22], may render the analysis of experimental data 
quite difficult. Using an extended self-similarity procedure (studying one moment against another one rather than as 
a function of the velocity difference) seem to provide better results for the inverse exponents, but since the cut-off 
effect is physical and not a measurement artifact, the exploitation of the data may lead to apparent contradictions 
|23]. 

It is to be noted that from duality, we may convert dual positive order moment to direct negative order moment. 
The latter may be very difficult to estimate experimentally because of the possible occurrence of arbitrarily small 
velocity differences over a given distance (or time using Taylor's hypothesis). Therefore, the duality relation may be 
exploited to obtain data which would be inaccessible otherwise. This procedure however relies on the applicability of 
this duality to experimental turbulence, which could not be tested directly if the direct moments cannot be computed. 
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FIG. 5: Scaling exponents obtained from simulations of the GOY shell model. The full curve gives the forward structure 
function exponents, r,. The dashed line connects the exponents 9q for the inverted structure functions. The black dot are the 
results obtained by applying the inversion formula Eq. 1221 and thus transferring the data for negative q's of to positive q's 
for Og. 

Nevertheless the present test of the analysis, using the shell model, constitutes an encouraging argument to proceed 
in this direction. 

Armed with these data sets we are ready to check the validity of the inversion formula (|22|l which by a simple 
inversion reads 

^riq) = 9-H-q) (28) 

Using the relation (|22|l directly by inserting T{q) for q G [—10; —1] we obtain by linear extrapolation the data for 
9q shown in Table II. The results are also shown as black dots in Fig. [S] Using the inverted relation (|28|) . we can 
compare the value of —Tq with the values of 9^^{—q), and these are also listed in Table II. Indeed, there is a very 
good correspondence between the values supporting the suggestion that the inversion formula is valid for the shell 
model turbulence data. This might be somewhat surprising as we do indeed perform a very different statistics. 

An inversion formula similar to the one derived in this paper has been proposed for the turbulence spectra 
and has been applied to multiaffine fields in 12^ . Recently, the inverse statistics has been applied to two-dimensional 
turbulence with the very interesting result that the inverse statistics of a smooth signal shows non-trivial behavior 
[26j | . Hastings has also derived a similar formula for Laplacian random walks in the very different context of diffusion 
limited aggregation . By using iterated conformal mappings Hastings obtain the exact multifractal spectra of the 
harmonic measure and derive the inversion formula for the f{a) spectrum. 

One might express a general worry that in the integrals for the statistical averages one integrates over the same 
singularity structure, both to obtain the standard and inverse (i.e. dual) structure functions. That is to say that 
dominating terms of the integrals (and thus the saddle point) come from the same singular structures even though one 
might argue that for turbulence, the velocity singularities are important for the forward structure functions whereas 
the laminar regions are important for the dual structure functions. However, as mentioned earlier, besides the intrinsic 
limitations due to the evolution of the cross-over length scale between laminar and inertial regimes with the moment 
order, when both tools are used to analyze the inertial regime, we are characterizing the same multifractal object for 
which the correspondence is expected to hold. Indeed the shell model used to validate the procedure is not deprived 
of such singularities (looking like shock waves in that case). 

VII. CONCLUSIONS 

The alternative approach to the standard multifractal spectrum and scaling exponents of different moment orders, 
which was proposed by interchanging the role of the physical quantity of interest and the length (or time) scale, has 
been examined in the case of a simple multifractal Cantor set. This example shows that the two spectra are simply 
related, and that the scaling exponents of the length moments can be related to the usual series of scaling exponents. 
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Based on this correspondence, we tested the apphcation of this ahernativc approach to turbulence using data 
obtained from sheU model calculations. A good agreement was found, thus suggesting that the above duality relation 
could be extended to more general cases. 
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TABLE I: Value of some r(g) and 9{q) exponents for moments of order q based on numerical simulations of the GOY shell 
model [iflL lid ITtI . These data are as presented in Ref. 0. For comparison, we present the corresponding series of exponents 
based on She and Levesque formula (index SL). To obtain the series of the inverse exponents 9sL{q) we have invoked the 
inversion formula 1221 . 
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TABLE IL Value of exponents obtained from simulations of the GOY model for negative moments in the interval q € [—12; —2]. 
Shown are the values of —T{q) which according to the inversion formula should be compared to 0~^{—q). Last row are values 
of 9{—T{q)) which should be compared to — g. 
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